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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 3254 


DETERMINATION OF FLAME TEMPERATURES FROM 20009 ТО 50009 K 
BY MICROWAVE ABSORPTION 


By Perry W. Kuhns 


SUMMARY 


Equations аге derived for the measurement of flame temperatures 
from the attenuation of a microwave beam by temperature-induced free 
electrons from chemical elements introduced in the flame. Corrections 
are developed to account for the electron distribution in the flame. 
Procedure for obtaining the temperature from absorption is outlined. 


The free electron collision frequency and an effective ionization 
potential for four alkali elements were determined experimentally. The 
qata were taken on а gas burner with a temperature range from 1900? to 
24009 K. The effective ionization potential of sodium agrees with the 
spectral-line limit value. From these experimental ionization potentials 
and the experimental collision frequency, ап accuracy of +600 К in the 
temperature was obtained. 


Fluctuating and average temperature data of a liquid propellant 
burner are presented primarily from 1.25-centimeter-wavelength microwave 
measurements in comparison with simultaneous two-color pyrometer and 
sound intensity measurements in the region from 22009 to 29009 к. 


INTRODUCTION 


Present experimental work on rockets and jet engines has shown the 
need for a method of measuring gas temperatures above 20000 K. Such а 
method to be practicable must be relatively insensitive to conditions 
outside the flame region of interest and to variations in gas flow. 
Because of the temperatures involved, such a method should not necessi- 
tate the presence of devices in the flame itself. Among such methods 
presently available are the sodium-line reversal, the "two-color" sys- 
tems, the infrared and the ultraviolet radiation methods, and the method 
making use of temperature-dependent attenuation of radiation in the 
microwave region. 
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The microwave attenuation method is dependent upon the absorption 
of electromagnetic waves by temperature-induced free electrons. The 
electron density, and thus the absorption power, is a function of the 
gas temperature. Іп most flames it is necessary to add an alkali salt 
of low ionization potential to supply sufficient free electrons. 


Belcher and Sugden (refs. 1 and 2) used microwave absorption to 
determine the electron density in & study of the ionization of alkali 
salts in a flame. They concluded that deviations of electron density 
from the values calculated from Sahats equation could be explained by 
reactions producing alkali hydroxide and oxide compounds, and the 
production of electron acceptors. Rudlin (ref. 3) presented data showing 
the variation of &ttenuation with flame temperature in the region of 
19009 to 25009 K and concluded that the attenuation was substantially 
independent of tbe halide component of the salt used. 


The present work is concerned with (1) the determination of the 
collision frequency of free electrons in & flame, (2) the estimation of 
the effect of nonequilibrium and flame chemistry upon the effective 
ionization potential of the alkali element, (3) the extension of the 
idealized infinite-sheet formula to the case of attenuation of a micro- 
wave beam by & finite cylindrical flame, (4) the estimation of sensitiv- 
ity of the attenuation measurement to variations in flame parameters 
such as temperature, flame size, gas flow, and salt concentration, and 
(5) the applicability of ап instrument for the measurement of average 
flame temperatures and of time-fluctuating flame temperatures. 


Microwave attenuation readings were taken on а blast burner for two 
wavelengths at temperatures between 19009 and 2400? K in conjunction with 
sodium-line reversal temperature measurements 2 and оп а liquid propellant 
burner at temperatures between 22009 and 2900“ K in conjunction with 
sodium-line reversal or two-color pyrometer measurements. Measurements 
of fleme-temperature fluctuations on the liquid propellant burner were 
taken with K-band microwaves in conjunction with two-color pyrometer and 


sound-intensity measurements. 
THEORY 


Unbounded Absorber 


Electromagnetic waves passing through an unbounded isotropic medium 
of free charged particles are attenuated exponentially: 


pepe ре (1) 


(А11 symbols used in the text are defined in appendix А.) 
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From classical electromagnetic theory (refs. 4 and 5 and appendix 
B), it is found that y, expressed in decibels per centimeter, may be 
closely approximated by: 


0.455 
reas (2) 
о + g€ 


where n is іп electrons per cubic centimeter and о and g are in 
reciprocal seconds. 


Margenau (ref. 6), using an energy-distribution equation instead of 
the mono-energic differential equation, obtains substantially- the same 
equation in the region of radiation frequency and electron density of 
interest as was obtained in this paper. 


In equation (2), g is a viscous damping factor due to the colli- 
sions between electrons and molecules. With the use of equation (2), 
the value of g can be determined experimentally by simultaneous meas- 
urements of the attenuation at two wavelengths 


(3) 


This method was used to determine g in this report. Classical theory 
sets g equal to the electron-molecule collision frequency (refs. 1, 2, 
and 4) which is found from the kinetic theory of gases (ref. 7) as: 


л 
g = pd” лет (4) 


Additional damping due to ion-electron collisions was calculated from 


the formula of Smerd and Westfold (ref. 8) and was found to be negligible. 


The electron density n ОҒ equation (2) can be determined from 
Saha's equation (refs. 7 and 9, and appendix C) assuming a single alkali 
salt in an electrically neutral flame: 


log Wd log f = 15.385 + 3 log T - =т= (5) 


2 
where п and N are in reciprocal cubic centimeters, а; is in electron 


volts, and T is in degrees Kelvin. А graph of n/N against log £/N 
is given in figure 1. In the region log f/N < -2, n/N= 4/f/N. 
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From equations (2) and (5) the variation of the attenuation coeffi- 
cient Y with temperature T may be obtained. This is demonstrated in 
figure 2 where the variation of y with T is shown for different con- 
centrations of sodium salt and for the K-band microwaves. 


Tt is seen from figure 1 that almost all atoms are ionized when 
log f/N >1. Above this value the change in attenuation with temperature 
becomes negligible. Тһе temperature at which f/N = 10 шау be described 
as the upper limit of usefulness of this method for measuring the tem- 
perature of unbounded flames. TIn addition, there is an instrument limi- 
tation in detecting microwave power which sets practical limits on the 
applicable electron density at 3x10l1/R < n< 5х1015/Ев, when R is in 
centimeters n is in reciprocal cubic centimeters. 


In the case of disequilibrium, Saha's equation may still be used if 
& changed ionization potential is assumed. Belcher and Sugden (refs. 1 
and 2) have demonstrated the disequilibrium effects and have proposed 
side reactions which will modify Saha's equation in a complex way. 
Assuming for simplicity that the only change taking place is in the ion- 
ization potential, then in the region Ғ/М << 1 the relation between 
two ionization potentials and the corresponding attenuation coefficients 
(from eqs. (2) ала (5)) is: i | 


2520( - ) 


The deviation of the experimentally obtained ionization potential from 
the spectral-line limit value is a measure of the effect of disequilib- 


rium on the ionization process. 


Finite Absorber 


In equation (2) the &bsorber is considered unbounded. Since the 
flame is of finite size, however, two effects will take place: (1) 
internal reflections due to the discontinuity at the flame boundary (ref. 
4), and (2) the effect due to the electron distribution in and the finite 
circular geometry of the attenuator. The first of these effects was 
checked theoretically by use of the equations in reference 4 and experi- 
mentally by measurements of the reflected beam intensities; the effect 
was found to be negligible. However, the effects of electron distribu- 
tion, a finite cylindrical absorber, and a finite beam cannot be 


neglected. 
Assuming a cylindrically symmetric attenuator traversed normally by 


a beam of square symmetry, the corrected attenuation coefficient Ae, ав 
a function of the parameters y and x/R, is defined by: 
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total transmitted pover 

атлар gim 10 deg y incident power (7) 
where A, is in decibels, x is one-half the mean width of the micro- 
wave beam, and R is the specific measured radius ag shown in figure 3. 


The corrected attenuation was calculated assuming the electron 
density to be a function of the flame radial distance. Shown in fig- 
ure 3(a) is the electron density distribution for the specific case of 
the blast burner determined from electrical probe and corrected thermo- 
couple probe measurements. Figure 3(b) shows a more general case assumed 
for a turbulent flame in thermal equilibrium. The microwave beam power 
pattern used is shown in figure 4, as determined by actual measurements 
for a beam emanating from a double concave lens and horn combination. 
This pattern was found to have almost the same shape in both the E and 
the H planes. The center line of the microwave beam was assumed to 
pass through the center of the cylinder and the direction of propagation 
was taken as normal to the axis of cylindrical symmetry. Because of the 
complexity of the problem, the integration was done numerically. The 
Plame width, beam width, and flame conductivity are fixed by the flame 
geometry and by. the given values of the parameters A and x/R. The 
quantity A is determined by assuming a homogeneous distribution of salt 
and a uniform temperature. The curves in figure 5 correct the attenua- 
tion for both the salt distribution and the flame geometry. The result- 
ing relation among A.s A, and x/R is shown in figure 5. 


For the finite absorber, the value of РИК for the upper Limit of 
usefulness for temperature measurement is more complex than the limit for 
an infinite absorber. A fair approximation taken from plotting corrected 


attenuation curves may be given as ?/М s 10 . 


Temperature Measurement Sensitivity 


With the insertion of values in the equations shown in appendix D, 
it was found that in the region used the sensitivity of the instrument 
to temperature was 4 to 15 times larger than the sensitivity to any other 
flame parameter such as flame radius, salt flow, flame velocity, or beam 
size- Since the coefficients B and C_ (appendix D) as shown in fig- 
ure 6 are larger for greater values of х/в, the final effect is to lower 
the temperature sensitivity as x/R increases until, in the region 
x/R = 1.8, the instrument is no longer a good temperature-sensing device. 


Because of the similarity іп slope and spacing of the correction 
curves in figure 5, the parameters B and C derived from figure 5(b) 
may also be employed to estimate the sensitivity of the temperature meas- 
urements on the blast burner flame. 
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| If a temperature change AT rather than the absolute temperature 
T is to be measured, an approximate expression for the uncertainty in 
such a measurement is given by: 


и ЕК + Білігі (8) 


More exact expressions for this uncertainty аге given іп appendix D. 


To simplify the reduction of data, an additional requirement for 
the measurement of fluctuating temperatures is that the change in inten- 
sity with temperature change be as linear as possible over the tempera- 
ture range. The extent of this linearity is demonstrated in figure 7 
where graphs of the attenuation as a function of temperature are plotted 
for five elements. 


The response of the system to time-fluctuating temperatures will 
depend upon the frequency and the velocity of propagation of these fluc- 
tuations. Assuming that the attenuation fluctuation moves at the gas- 
stream velocity u and varies sinusoidally with time at frequency v, 
the frequency v2 at which the indication of amplitude of the alternat- 


ing components of power will be reduced 1 decibel is given by: 


BITS 


EXPERIMENTAL METHODS 
Procedure for Temperature Measurement 


Attenuation readings were taken of the flame with and without a 
known amount of salt added. ТР the attenuation without added salt could 
not be neglected (which happened only with the alcohol-oxygen flame), 
the attenuation with added salt was corrected by means of equation (C2). 
The temperature was found graphically from curves of attenuation against 
temperature obtained from equations (5) and (2) and figures 1 and 5. An 
example showing the steps used may be found in appendix E. · 


Instrumentation 


A schematic diagram of the measuring system is shown in figure 8. 
The microwave generator for the X-band measurements was a U.S. Navy 
Т5-120 test set containing a 723-A/B Klystron oscillating at 9650470 
megacycles. The microwave generator for the K-band measurements was 
@ power supply of NACA construction with a 2K 33 Klystron oscillating at 
27,600+100 megacycles. The reflector voltages of both Klystrons were 


modulated. 
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Measurements of the beam intensity patterns showed that in order to 
make the system applicable to flames of the size used, it would be nec- 
essary to narrow the beam by means of double-concave metal converging 
lenses. The lenses used were of aluminum sheet in Lucite holders, de- 
signed from information in reference 10. The values of х obtained were 
1-8 centimeters for the K-band and 3.2 centimeters for the X-band. The 
distance between lenses was approximately 40 centimeters, with the flame 
midway between lenses. Attenuation of the beam by ionized flames did 
not change the voltage standing wave ratio by any detectable amount. 


Silicon diode crystals were used as square-law detectors. А 1N26 
crystal was used for the K-band (1.25-cm wavelength) and a 1N23 crystal 
for the X-band (3-cm wavelength) measurements. Two methods were used at 
various times for measuring the crystal outputs: a tuned voltmeter 
reading in decibels (with 1 kilocycle modulation on the Klystron) for 
measuring mean temperatures; and an amplifier, band-pass filter, and 
cathode ray oscilloscope system (with 20-kilocycle modulation on the 
Klystron) for recording varying flame temperatures. 


Blast-Burner Tests 


A blast burner which burned propane-air-oxygen mixtures was used 
for the experiments. By varying the fuel-oxidant ratio. the temperature 
could be varied from 1900° to 2500° K. The burner was mounted on a 
large optical bench in such a way that it could be moved to intercept 
either the K-band or the X-band beams. 


The alkali elements were introduced as water solutions of halide 
salts in a fine spray at the base of the flame. Readings of attenuation 
by the flame were made with and without the spray. The atomic flow rate 
was found by measuring the weight of the solution used over a given 
period of time, the flow rate being kept constant. The flame volume flow 
rate was found from integrations of total-pressure measurements taken 
across a number of flames. The average atomic density was obtained by 
dividing the atomic flow rate by the flame volume flow rate. As the 
salt was not uniformly distributed throughout the flame, a method was 
devised by which the magnitude and the shape of the electron density 
curve in the region below 21000 K could be estimated. An ionization - 
probe with a separation of 2 centimeters between poles was moved through 
the flame and both current and voltage were measured; from reference 11, 
n © i/U,V. Because of the unknown values of electron mobility and field 
strength only the relative electron densities could be determined. At 
the same time, thermocouple temperature measurements were taken to deter- 
mine the relative temperature profile of the flame. То correct for radi- 
&tion effects, the thermocouple readings were raised until the peak 
reading equaled the sodium-line reversal temperature. From these data a 
curve of the relative atomic density could be obtained. By integrating 
and normalizing this curve over the flame cross section, a value was 
obtained for the ratio of. the peak electron density to the average 
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electron density computed from the salt flow rate for the assumption о? 
uniform temperature and salt distribution. 


The width of the flame was measured in three different ways: (1) 
by means of a platinum probe connected to a commercial electronic switch 
which closed when the resistance between the probe ends was less than 
2 megobms; (2) from total-pressure measurements with a probe; and (3) 
by means of a salt-coated pyrex tube, the salt being melted from the 
tube when placed across the flame. 


To ascertain the temperature of the flame by the sodium-line re- 
versal method (ref. 12), the sodium was introduced as NaNOz powder at 


the flame base. 


The blast burner was used to determine the electron collision fre- 
quency by use of equation (3) and the ionization potentials by use of 
equation (6) and to check the correction technique. 


Liquid-Propellant Burner 


A photograph of the apparatus in the burner cell is shown in fig- 
ure 9. 


The propellants (alcohol-oxygen, and heptane-oxygen that contained a 
small amount of turpentine) were ignited in a tube 5.1 centimeters in 
diameter which was open at one end. The temperature range of 2200° to 
2900° K was obtained by varying the oxidant-fuel mass flow ratio in 
discrete steps of 1.2, 1.6, 2.0, and 2.4. The propellant flow was held 
constant at 0.2 pound per second. The microwave beam was "focused" in 
the exhaust, 8 to 18 centimeters below the tube end. Because of space 
limitations, it was possible to use only one wavelength at a time. 

Both wavelengths were used in measurements on the &lcohol-oxygen sys- 
tem, while only K-band measurements were taken on the heptane-oxygen 


system. 


The liquid-propellant burner was used to establish the feasibility 
of using the microwave apparatus on such types of flames and to obtain 
data on the rapid fluctuation of temperatures in flames. 


The salt used, Nal, was introduced directly into the fuel. The 
salt was easily soluble in the alcohol, but not in heptane. Thus, it 
was necessary first to dissolve the NaI in acetone and then to suspend 
this solution in the heptane. Other alkali salts were tried with 
unsatisfactory results as to mixing. Since from prior experiments an 
unknown amount of sodium salt contaminated the alcohol-oxygen flame, 


equation (C2) was used in determining the temperature. 


nee ee Oe ----- 
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Temperature measurements were also made in the upper part of the 
tube for an oxygen-alcohol flame by sodium-line reversal (ref. 13) and 
at the level of the microwave beam for the heptane-oxygen flame by the 
two-color pyrometer system (ref. 14). 


EXPERIMENTAL RESULIS 
Blast Burner 


Determination of correction for electron distribution. - Typical 
data taken with an ionization probe and with a thermocouple probe are 
shown in figure 10. From averages of such curves the shape of the 
electron-density curve (fig. 3(a)) was determined. The absolute magni- 
tude of the curve was determined by equating the area integral of the 
normalized electron density curve to the average electron density cal- 
culated from the measured gas flow, salt flow, and sodium-line reversal 
temperature. The ratio of the maximum electron density to the electron 
density, assuming a homogeneous distribution of salt, was found to be 
0.55. From this curve, the geometry correction for the specific case 
of the blast burner flame has been determined as shown in figure 5(a). 


A partial experimental verification of equation (8) for absorption 
by a finite-size flame was found by means of measurements in the X-band 
region with no lenses, one lens, and two lenses. The ratios of experi- 
mentally measured attenuations for these three conditions were 
0.64:0.77:1.0, respectively, while calculations by means of figure 5(a) 
predicted ratios of 0.40:0.80:1.0. The discrepancy in the results 
when no lenses were used is due to the inability of the correction to 
compensate fully for such a large ratio of x/R: 


Determination of collision frequency. - From simultaneous attenua- 
tion readings on the blast burner in the region T = 22009 K, the ratio 
of y for the X-band microwaves to y for the K-band microwaves, after 
correction for beam size, was found to be equal to 1.65. Insertion of 
this value in equation (3) along with known values of the microwave fre- 
quencies gave a value of д=1.9х1011 sec-l. This value may be compared 
with the kinetic theory value of g#0.8x1011 sec-l (eq. 3), the one 
determined from the electron mobility (eq. (В14)) as used by Wilson 
(refs. 11 and 15) in flame ionization studies of 6%5.9Х1011 sec-l, and 
the one determined experimentally from microwave attenuations by Belcher 
and Sugden (refs. 1 and 2) of gc«0.9»30ll sec-l. Іп the region 
0.8=1011 < g < 2х1011, K-band attenuation readings are fairly insensi- 
tive to the collision frequency g, while the X-band attenuation in this 
same region is approximatly proportional to 1/ Ze 


Іп extending the value of collision frequency to other temperature 
regions beyond that in which measurements were made, it was assumed that 


g e p/A/T ав in equation (4). 
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Determination of effective ienization potential. - The total- 
pressure profiles of 27 different flames were measured. From these pro- " 
files and the temperature profile (fig. 10), the integrated volume flow 
rate was obtained as a function of the sodium-line reversal temperature 
(fig. 11(а)). The value of 2R was found from flame-width measurement 
by three methods and is shown in figure 11(b) as 5.7 centimeters. 


The theoretical attenuation coefficient at 20009 K for the accepted 
ionization potentials (ref. 16) was calculated. From average values of 
the K-band attenuation at 20009 K, the effective ionization potentials 
were found from equation (6) and are shown for four elements (cesium, 
rubidium, potassium, and sodium) in table I. Also shown are the values 
found by Wilson (ref. 16) and Belcher and Sugden (ref. 2) using the g 
employed in their writings and using the g of this report. The effect 
of changing the g is to raise Wilson's results by about 0.4 ev and to 
lower Belcher and Sugden's value by about 0.2 ev. 


5109 _ 


It can be seen from table I that sodium seems unaffected by flame 
chemistry or disequilibrium and that the stronger the binding of the 


electron the smaller are the effects. 


From the experimentally obtained ionization potentials for cesium, 
rubidium, and potassium and the accepted value for sodium, the micro- 
wave temperatures were obtained by the procedure outlined in appendix 
E. The resultant temperatures are plotted against the sodium-line 
temperatures in figure 12. 


From figures 11(а) and 11(b) the probable errors in atom density 
and flame width are AN/N=+0-10 and AR/R=i0.15. The uncertainty in 
the sodium-line reversal temperature was taken ав AT/T=+0.01. When 
these values and the values of the parameters qj, B, and C are 


inserted into equation (D6), the average deviation expected between 
sodium-line reversal temperatures and microwave temperatures should 
range from about 559 K at 19000 К to about +70° K at 24000 K. Of the 
K-band and X-band points, 81 and 67 percent, respectively, in fig- 

ure 12 fall within these limits. However, there is a small systematic 
error in the X-band tests. It is suggested that this error is probably 
due either to failure of the theory to correct the attenuation completely 
for flame geometry in the higher attenuation regions or to the strong: 
dependence of the X-band attenuation on the electron collision frequency. 
Thus some of the advantage that an X-band system enjoys due to larger 
and more available components is offset by the larger beam size and the 
collision frequency dependence. 


As is seen from figure 12, there is no apparent dependence of the 
ionization potential upon the oxygen-nitrogen ratio in the flame. It 
can be concluded that sodium can be used in a hydrocarbon flame with 
the accepted value of the ionization potential in Saha's equation. 
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Liquid Propellant Burner 


Mean temperature measurements. - In the case of the liquid propel- 
lant burner, fewer individual measurements were made than for the other 
burners. The flame width measured from photographs was 7.6 centimeters. 
The flame velocity calculated from high-speed photographs of a screaming 
condition and from the frequency response (eq. (9)) was 3*10* centi- 
meters per second, giving a flow rate of 1.36106 cubic centimeters per 
second. If an average molecular weight of the exhaust gas of l5 is 
assumed, з flow rate of 1.2x106 cubic centimeters per second is found 
from the propellant flow rate. A value of 1.5Х106 cubic centimeters per 
second was used in the computations. The amount of salt placed in the 
fuel varied from 0.4 to 1.0 grams in 4000 cubic centimeters of fuel. 
Atomic densities in the exnaust varied from 1.0 to 1.7x1015 atoms per 


- cubic centimeter. The temperature measurements on the oxygen-alcohol 


burner shown in table II were found by averaging the results of micro- 
wave voltmeter readings of numerous runs for given fuel-oxygen mass-flow 
ratios. The average difference between the two results is 60° K. Use 

of the classical collision frequency of g= 0.8х1011 raises the differ- 
ence between the two results to about 200° K, well outside the experi- 
mental error. On the oxygen-heptane burner, averages over short, 
randomly-spaced intervals were taken by use of the cathode-ray oscillo- 
scope apparatus. The results for K-band microwaves are compared with 
average "two-color" measurements in figure 15. 


Measurement of temperature fluctuations. - When the crystal output 
was connected to the plates of a cathode-ray oscilloscope, the rapid 
fluctuations in the attenuation could be observed. Use of equation (9) 
shows that X-band microwaves would be unsatisfactory in measuring fluc- 
tuating temperatures in this flame; therefore only K-band data were 
taken. A test was made to decide whether the fluctuations were primarily 
due to temperature changes or to salt density changes. The test was 
based on the fact that, with a low concentration of cesium salt in the 
flame at the temperature of the liquid-propellant burner flame, the salt 
would be almost totally ionized and any large variation in the attenua- 
tion would be attributable to flow variations. This experiment was 
performed and yielded only small attenuation variations at those fuel 
mixture settings at which injection of sodium salt yielded violent atten- 
uation fluctuations. These results indicated that the observed attenua- 
tion fluctuations were due almost entirely to temperature fluctuations 
and not to changes in the gas flow rate. 


Figure 14 shows measurements on the heptane-oxygen burner of temper- 
atures obtained by K-band microwave attenuation measurements, and by the 
two-color pyrometer. Also shown are simultaneous sound-intensity meas- 
urements taken with a microphone 15 centimeters away from the burner. As 
the response of the microphone for the conditions used is not definitely 
known, the microphone data are included only to show that screaming 
conditions existed during the measurement. 
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If both the mean temperature and the measured change in attenuation 
are estimated to be accurate within 2 percent, the error in the change 
in temperature is 4.5 percent. Figures 14(a) aad 14(b) represent condi- 
tions of screaming at 1000 and 3000 cycles per second, respectively. 

The two-color instrument shows higher frequency fluctuations than the 

microwave instrument because the latter averages over a volume 100 times 
as great as that covered by the former. Two longer measurement periods 
are shown in figure 15 for the sound and K-band microwaves. Additional 


data may be found in reference 14: 


In the figures the flow rate and the flame radius were assumed 
constant for the microwave data, and no correction was made for the fre- 
quency response to attenuation fluctuations. Taking the effect of these 
into account, it is estimated that for 1000 cycles per second screaming 
the temperature peaks are 20° K low, the valleys 20° K high, and no 
change occurs in the rest of the curve. 


Effects of Addition of Salt Solutions on Flames 


The addition of the water-salt solution to the blast flame was 
Observed to change the region of flame stability and to change the flame 
sound. The effects could nearly be duplicated by using a pure water 


spray. 


When alcohol and oxygen were used in the liquid-propellant burner, 
the addition of salt had no noticeable effect on the flame except to 
color it yellow. In the case of heptane-oxygen, the addition of acetone 
in which sodium iodide was dissolved was found to һауе a decided effect. 
With no acetone, steady screaming was observed only at an oxygen-fuel 
mass-flow ratio O/F = 1.6; with 3 percent acetone in the fuel, scream- 
ing was observed at О/Е = 1.6 and 2.0. With 8 percent acetone, the 
burner screamed at O/F = 1.2, 1.6, 2.0, and 2.4. Тһе amount of acetone 
used in the temperature-determination experiments was always 3 percent 


or less. 
CONCLUSIONS 


The following conclusions шау be drawn from the experimental and 
theoretical investigation of flame temperatures by microwave absorption: 


1. Theoretical computations indicate that microwave attenuations 
can be used as a measure of both mean and rapidly fluctuating tempera- 


tures, with little dependence upon gas flow. Experimental studies verify 


these computations within the error of the measurements. Over the range 
used, the theoretical sensitivity of the attenuation to temperature was 
estimated to be 5 to 15 times larger than the sensitivity to any other 


flame parameter. 
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2. The experimentally obtained collision frequency for electrons 
was found to be approximately 1.9х1014 reciprocal seconds at a tempera- 
ture of 22009 K. 


3. The ionization of sodium seems unaffected by flame chemistry, 
and the accepted value of the ionization potential may be used. The 
elements cesium, rubidium, and potassium may be affected by flame 
chemistry in such a way as to necessitate a calibration by adjusting 
the values of the ionization potentials. The ionization potentials seem 
to be independent of the nitrogen-oxygen ratio of the flame. 


4. The attenuation formula for an unbounded sheet has been extended 
to include finite flames. Suitable corrections can be made for the size 
of the flame if the flame diameter is 60 percent or more of the beam 
width. Smaller flames introduce serious errors. 


5. Microwave temperatures can be read from oscilloscope records by 
use of a reference scale. Over a temperature span of 600° K (the mid- 
point depending on the salt used) the variation of intensity with tem- 
perature is almost linear. Over this linear range the temperature will 
be accurate to within about 2609 K and changes in temperature can be 
measured within +50 K or 4.5 percent of the temperature difference, 
whichever is larger. The range of the instrument may be increased by 
the use of & meter with linear decibel scale. 


6. The use of K-band radiation was advantageous because of the 
smaller beam size and the insensitivity of the attenuation to electron 
collision frequencies. 


7. The effect of the addition of the required quantity of salt to 
the flame has no effect on the combustion properties. However, if addi- 
tional material is added as a solvent for the salt, the burning may be 
effected drastically. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 26, 1954 
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SYMBOLS 
The following symbols are used in this report: 
attenuation, A = 2Ry = zy 
corrected attenuation, A, = ZRT o 
B,C parameters 
c velocity of light in vacuum, 2.998x109 m/sec 


d molecular diameter 


E electric field strength along the y axis 
e electronic charge, 1.602х10-19 coulomb 
f function in Saha's equation 


viscous damping coefficient (collision frequency) 


g 

H magnetic field strength 

І intensity of the microwave beam 
i 


ionization probe current 


J current density 

j Na 

кю Boltzmann's constant, 8.628Х1079 ev/°K ог 1.580х10716 erg/K 
m electron mass, 9.107х10-51 kg 

N number of atoms (per unit volume) furnishing free electrons 

n number of free electrons per unit volume 


ny number of positive ions per unit volume 


т 


microwave power 


p gas static pressure 
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х,у,2 


ionization potential 

specific measured flame radius (fig. 2) 

radial distance vector 

absolute temperature 

time 

electron mobility 

flame velocity 

ionization probe potential 

salt mass flow rate 

distance vectors 

one-half mean width of the microwave beam 

phase constant 

attenuation constant 

power attenuation coefficient for an infinite flame 
power attenuation coefficient for finite flame, db/cm 
electrical permittivity 

permittivity of free PRSE. 8.85x10-1¢ farad/meter 
propagation constant, x = o + if 

magnetic permeability 

permeability of free space, 1.257х106 henry/meter 
cyclic frequency of fluctuation 

cyclic frequency at which the measurement is down 1 db 
electrical conductivity 


angular frequency of radiation 
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Subscripts: 
O initial conditions 


1,2,5 used to denote separate conditions appearing in the same 


equation 


— —À — — 
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DEVELOPMENT OF ATTENUATION EQUATIONS 
The following equations are developed in the same sequence as those 
in reference 4. А11 units of electrical measurement are rationalized 
units. 


The differential equation for a plane wave propagating in the z- 
direction in an unbounded medium is: 


2 2 
OTE ae ОЁ sin OR sg 


3:2 ot? 9% еч 


Тһе equation Рог H is identical. 


ІҒ the electromagnetic wave is harmonic in time propagating in the 
positive z-direction, then the solution of (B1) is given by: 


E = Eg exp (-jot + jaz - 82) `_ (B2) 
where 
х2 = pea? + juco = (a + jp)? (85) 


and a is the phase constant and B is the amplitude attenuation 
constant. 


The equation for the displacement of а free electron in а polarized 
electric field. is: 


Š | 
m 25 + mg SY = e Eo exp (-jot;) (B4) 


where y is the displacement of the electron from equilibrium and g 
is a viscous damping factor. Тһе steady-state solution of (B4) is then: 


y = — 1 ej Ес exp ( - just ) (B5) 
g - jo om 


If there are n free electrons per unit volume, the maximum current 
density J is: 


2 
z [пе ДЕ 
I= пе H (пе) (в6) 
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since J = ОЕ, it follows that 


glo, 


n 


б (В7) 


 &- jo 


Assuming both 6 = So and р = ро, insertion of the value of G 
from (B7) into equation (B3) yields B as: 


B56 /2 
poii.) 044 в a 5 (88) 
C 2 фе (1 Ж Ф) 
where 
2 
p = Дес 3) (во) 
ME, (2 + E? 
For Ф<<1, by means of the binomial expansion, 
p = EÈ (1 + z) (B10) 
26 2 


In the wavelength and concentration range of interest, B сап be 
further approximated by 


AE = --- amu. M" (в11) 
Since 
E осет?“ 
or 


P= Po е SPS = Po e f? 


insertion of numerical values for e, c, m, and По in the proper units 


yields 


үс 0.455 (в12) 
"es 
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where ү is in decibels рег centimeter, n is іп electrons per cubic 
centimeter, and g is in reciprocal seconds. Equation (B12) may be used 
up to the region y “5 for K-band microwaves and ү * 10 decibels per 
centimeter for  X-band microwaves. 


The d-c mobility is defined in reference 7 as 


02 = A . 5 (B13) 
Substituting from equation (B6) and setting 6>>ш 
e 
Ue “ cu (B14) 


— — — 
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SAHA'S EQUATION 


Saha's equation (refs. 7 and 9) for alkali elements is given by: 


E dr EM 15/2 exp (хт) (с1) 


N- n \ n2 KT 


For a single salt in an electrically neutral flame n, = п and equa- 
tion (C1) may be expressed as equation (5). 


Of interest is the problem of attenuation measurements of a flame 
which is contaminated by an unknown quantity of salt. Іп this case, two 
attenuations may be measured, y, and Yo, the respective attenuations 


before and after the addition of a known amount of salt. The following 
equations may be set up: 


for contaminating salt in the flame 


n11 


жые =: 
Ny - пуу 1 


for contaminating salt plus known salt 
(nig + nee) ше, 
N. - n ШЕ; 
qum QE 
(012 + ngg)(mig + n22 - ча) 
No - 122 
for known salt only in the flame 
2 
= Ов 
No - nog 2 


where the first subscript refers to the salt and the second to the at- 
tenuation reading. ü 


fo 


The + 0202 © Vo 


nj @% Ty 
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| Defining 
0.455n528 
g^ + оё 
and defining 
fo > о 
ry + 13 - 1g = 0 (са) 
if f4 = fo; then 
a ы _ 0.4558. | 
Үр - Тү 
2 ос. y2 aw ye [212 "4 
үү deye ys 13 ( та ) (C4) 


mn ie „ „м U U = 
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APPENDIX D 


SENSITIVITY OF TEMPERATURE MEASUREMENT 


The following approximate formulas serve to indicate the sensi- 
tivity of the microwave temperature method to changes in various para- 
meters and to aid in estimating possible errors. Equations (6) and (2) 
be cambined to give the dependency of y on T over the region 
P/N < 1. 


Y © Ба т1/4 on(- ЕЈ Q^ + g^ (01) 
This dependency may be approximately represented by the proportionality 
Y о A exp (- == (D2) 

also 
Y = Y. Tre (05) 


From equations (Dl) and (D2), 


тю (2x 5а) | (D4) 


ə(r/T. a(r/r, 
ss Z а. s 


where А = 2Кү. Combining these equations yields the following relation 
among the variations in the experimental factors which enter into a 
temperature determination: 


and 


— m а= == 


әт „кт 2 97 SN, „В- Сок, 20 _ 5x (гв) 
ї-вү Ñ 1-BR 1-387 


where the parameters 
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С = Ә(т/те) | хд 
Ə(x/R) Т/Т 
are obtained from figure 5(b) and are shown in figure 6 as functions of 
A and x/R. ` Since 


and 


city u 


(D7) 


In the range of usefulness of a single salt, q /ET will vary from 20 


to 50. By inserting in equation (D7) the extreme values of B and С 
from figure 6, it is found that the range of the ratio of the sensitivity 
of the attenuation to temperature to the sensitivity of the attenuation 
to any other flame parameter is: 


Ratio of sensitivity 


From these values it can be said that in the range of attenuations where 
the instrument was used the attenuation was 5 to 15 times more sensitive 
to temperature change than to any other flame parameter. Тһе sensitivity 
to flame radius becomes larger with increasing ХИВ. From these values 
it can be seen that serious errors will enter when X/R is greater than 
1.4 and that above Х/Б = 1.8 the instrument is unsuited for the de- 
termination of temperature. 


If the temperature change AT from a mean temperature Т is to 
be measured by the deviation of the attenuation AM, from the mean 


attenuation A,, then the sensitivity of the change using equation (D6) 
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and assuming R, x, and М constant is given by 


s(at) _ ^%һ) жт 8А, 
ga cp (28) 
с T A. 


where T/T is to be substituted from (D6). 


The usual method of determining fluctuations in temperature is by 
use of previously prepared theoretical curves of intensity ratio against 
temperature for different values of N. If the mean temperature is 
determined from these curves and the measured microwave attenuation, 
then by use of equations (D7) and (D8) the probable error in the tem- 
perature fluctuations is given by: 


-- 16 
Bess „ Св) (E amy = ZO М 
AA, 2kT K. 
2 2 2 2 
2C SR Ow Su 20 Әх 
ER E Ы rx «RI TTBS (D9) 


where the 65's denote probable errors. 


If; however, an independent method of measuring the mean tempera- 
ture is used then equations (D6) and (D8) yield 


„(ој ИЕЛІ ЕЗЕЗЕ: 


ae s. + а - B) |, + [a s|: 


l-B'x 


Since in most cases the mean attenuation is taken from a large 
number of readings, the terms containing A, may be neglected. 


— --- —— —— ee —— . — — —_ 
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EXAMPLE FOR OBTAINING THE TEMPERATURE FROM THE ATTENUATION 


lf sodium salt and K-band attenuation are used, the experimental 
attenuation due to the added salt A. may be changed to A by means 


of figure 5 and the known value of x/R; A is divided by the flame 
diameter ZR to obtain y, апа Т is found from figure 2 for the 
known value of М. 


If another salt and/or another microwave frequency is used, the 
following example gives the procedure by which y is found as a 
function of T as in figure 2. 


L. The following conditions are assumed: 


Flame volume flow rate, сс/вес.................. 108 
Salt mass flow rate, gm/sec .................?2.2х075 
G) (K-band microwaves) оо оо е òo ө е e ө е 9 % е ее @ ө ds 73x10 


s Lal EU ЕЕ ЕВ 8.9жі07%///т 


Salt used e ө ° ° % ° е ° е е Ф Ф е е . Lj Ф е ° е * Ф . e е ° е 
q; 3 ev . . Ф e . . . . . . L4 ° е ° e ° ° е е е a Ф e ° ° е . ° е Js ғы 


Mole cular wel ght Ф ° ° ° 6 ° ә е е . * e ° е е ° + ° ° ° ° ө а ә 155 ° 9 


2. Find: ү at 28009 K 


E 2 
5. - 2:290 76.02 X0 XL. 1015 atoms of та ее 
106х155.9 
4. log f = 15.385 + 2 log 2800 - MN = 10.90 (eq. (5)) 


log f/N = -2.10, n/N = 0.089 (fig. 1) 


5. rt ‚в (eq. (2) 


15.Оха/М = 1.16, db/cm 


By repeating steps 4 and 5 for different temperatures, a graph similar to 
figure 2 пау be plotted. For known values of x and R, a graph of 
the attentuation A. against the temperature may be plotted with the 


use of figure 5. 


- —— n - — — 
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TABLE I. - COMPARISON OF IONIZATION POTENTIALS OF ALKALI METALS 


Damping Tonization potential, ev Estimated 


= ev 


Accepted, 
reference 16 
Present data 
Reference 15 
Reference 15 
Reference 2 
Reference 2 


TABLE ІІ. - MEAN MICROWAVE TEMPERATURES 
OF OXYGEN-ALCOHOL BURNER EXHAUST USING 


K-BAND AND X-BAND MICROWAVES 


Oxidant-fuel K-Band X-Band 
mass-flow ratio | CK OK 
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5 -2 -1 О 1 
log f/N (increasing temperature ——) 


-4 


ic (53) Ratio n/N ав a function of log f/N, by means of Saha's equation 
eq. . 
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gure 2. - Attenuation coefficient y ав а function of the temperature 
for K-band microwaves. Sodium salt; experimental value of g; different 


values of N. ү = 1.3x10-12N(n/N). 
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(a) Specific case of blast burner taken (b) General case, assumed for liquid 
fron measurements of temperature and propellant burner. 
ionization ourrent in flame, 


Figure 5. - Electron density ав a function of flame radial distance. Used in cal- 
culation of Plame geometry correction. 
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Figure 4. - Measured microwave beam intensity pattern as a function of 
lateral distance x from beam center line. Used in calculation of 
flame geometry correction. 
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(a) Flame geometry for specific case of blast burner (see fig. 5(а)). 


as a function of attenuation A 


- Corrected attenuation A, 


for values of the parameter X/R. 


Figure 5. 
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(b) Flame geometry for general case (see fig. 3(b)). 


as a function of 


Ao 


attenuation A for values of the parameter x/R. 


Figure 5. - Concluded. Corrected attenuation 
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Attenuation, A = 


- Coefficients B and C of the sensitivity equations 


against attenuation A with Y/R as parameter (fig. 5(5)). 


Figure 6. 
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Figure 7. - Theoretical microwave power attenuation for K-band radiation 


1015 atoms per cubic 


2R = 7.6 centimeters; М 


centimeter; g = 1.9104, / 2200; 44, accepted values. 


for five elements. 
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Figure 8. - Schematic diagram of the temperature-measuring apparatus. 
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Figure 9. - K-band microwave equipment in the liquid propellant burner test cell. 
(electronic equipment in controlroon.) 
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(b) Flame width against temperature by three methods of 
width measurement (2R = 5.7 сш). 


Figure 11. - Experimentally obtained flame parameters of 
blast burner against sodium line reversal temperatures. 
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(a) K-band microwaves, obtained for experimental ionization 
potentials for cesium, rubidium, and potassium. 


Figure 12. - Microwave temperature against sodium line reversal | 
temperature for blast burner. š 
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(b) K-band microwaves, obtained for accepted value of 
ionization potential of sodium (5.12 ev). 


Figure 12. - Continued. Microwave temperature against 
sodium line reversal temperature for blast burner. 
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(е) X-band microwaves, obtained for experimental ionization 
potentials for cesium, rubidium, and potassium. 


Figure 12. - Concluded. Microwave temperature against sodium line 
reversal temperature for blast burner. 
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Figure 15. - Comparison of mean temperatures found by two methods of 


measurement on oxygen-heptane burner for three oxygen-fuel mass- 


flow ratios. 
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(a) Screaming at approximately 1000 cycles per sec 


Figure 14. 


- Typical temperature variations in liquid propellant burner. 
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- Concluded. Typical temperature variations in liquid propellant burner. 
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